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Based on the Schrodinger equation, exact expressions for the non-relativistic 
particle energy in the local external field and the external field potential are 
derived as inhomogeneous density functionals. On this basis, it is shown that, 
when considering more than two noninteracting electrons, the energy of such a 
system cannot be an inhomogeneous density functional. The result is extended for 
the system of interacting electrons. This means that the Hohenberg-Kohn lemma 
which assert that in the ground state to each inhomogeneous density corresponds 
only one potential of the external field cannot be a justification of the existence of 
the universal density functional in the general case. At the same time, statements 
of the density functional theory remain valid when considering any number of 
noninteracting ground-state bosons due to the Bose condensation effect. 
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According to the Hohenberg-Kohn lemma [1,2], the same inhomogeneous density n(r) 
cannot correspond to two different local potentials t>i(r) and ^(r) of the external field in the 
ground state of the non-relativistic system of electrons (except the case Ui(r)— t> 2 (r) = const). 
Thus, the inhomogeneous density n(r) of the ground-state non-relativistic electron system 
uniquely defines the potential v(r) (to within an additive constant). In the case of ground 
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state degeneration, the lemma relates to the density n(r) of any ground state. It is generally 
accepted [1] that from the Hohenberg-Kohn lemma the external field potential v(r) is an 
inhomogeneous density functional, 

n(r) = n[v (r)] — > v(r) = v[n(r)] + const. (1) 

Existence of the functional (the square brackets in (1)) means that there is a universal rule 
according to which the external field potential v(r) can be found by the known inhomoge- 
neous density n(r) corresponding to the ground state of a system under consideration. This 
means that in principle there exists (although cannot be explicitly found or indicated) the 
rule of determining the function v(r) by the known function n(r) whose structure is inde- 
pendent of the explicit form of functions v(r) and n(r). We emphasize, that this very strong 
statement has not yet been called into question. In fact, in the general case, there is no one- 
to-one correspondence between the Hohenberg-Kohn lemma and relation (1) on functional 
dependence. It only follows from the Hohenberg-Kohn lemma that quite a definite external 
field v(r) can be put in correspondence with each function n(r) (within a constant factor). 
However, this does not mean that such a correspondence is established by the unified rule 
v(r) = v[n(r)] universal for any external field [3]. 

In other words, each external field determines the unique density (it is clear, e.g., on the 
basis of the universal rules of the perturbation theory), and each density determines the 
unique external field on the basis of the Hohenberg-Kohn lemma. However, the rule for 
the last correspondence can be non-universal. This rule in general depends on the concrete 
form of the density. Opportunity of such non-universality of the rule of correspondence for 
the "inverse" functional v(r) = v[n(r)] was not considered in [1]. Existence of this non- 
universality violates the Hohenberg-Kohn theorem, although the Hohenberg-Kohn lemma is 
undoubtedly correct. Below we show that the universal density functional does not exist in 
the case of the number of fermions more than two. 

For this reason, let us pay attention to the fact that the inhomogeneous density n(r) 
is the functional n[v (r)] by definition n(r) = (^o I x I /+ (r) ^(r) | \J/ ) and the functional 
dependence of the wave function \l/o of the ground state on the external field v(r), — 
ty [v(r)]. Here ^ + (r) and \&(r) are the field creation and annihilation operators. Thus, 
the correspondence rule between the density n(r) and field v(r) is established based on 
the solution of the corresponding Schrodinger equation for the wave function \I/ in a given 
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external field v(r) and the inhomogeneous density n(r) determination. Then, it only follows 
from the lemma proved by Hohenberg and Kohn [1] that the functional n[t>(r)] is unique 
(taking into account the condition n[w(r)] = n[v(r) + const]). Here it is clear that the 
functional n[i>(r)] is essentially nonlinear in the external field v(r). This means that two 
possibilities are admissible without violating the Hohenberg-Kohn lemma: (i) the inverse 
problem on the determination of the dependence of v (r) on n(r) has individual solutions 
to each pair of functions n(r) and v(r) (or for certain types (classes) of pairs of functions 
n(r) and v(r)); (ii) the inverse problem has the universal solution v(r) = v[n(r)]. As noted 
above, this dilemma is not usually considered, and it is assumed that there is the universal 
solution v [n(r)] valid for any external field and any number of particles, i.e., the possibility 
(ii) is always realized [3]. 

In fact, as follows from the above, the Hohenberg-Kohn lemma is insufficient for the 
statement about the existence of the universal solution t>[n(r)]. However, it seems impossible 
to disprove the statement about the universality in the general form. In this regard, we use 
the "proof to the contrary". Let us assume that the functional v(r) = v[n(r)] exists and 
analyze consequences of this statement. By the example of noninteracting fermions, let us 
show that such an assumption leads to contradiction. 

If we accept the validity of statement (1), then the energy E the ground-state system 
of N interacting electrons with the Hamiltonian H in the external field with potential v(r), 
which is characterized by the wave function ty , can be written as 



E = <*o|# |*o> = E (N, [*oXr)D = E [n(r),v(r)}, N = / n(r)ci 3 r. (2) 



Here it is taken into account that (^ |^o) = 1 an d \l/o[^(r)] = ^o[v(r) + const]. In turn, 
it immediately follows from (2) that the quantity F[n(r)] = (^o\T + U\^o) which defines 
the system ground state energy 



is the functional of only the density n(r) ("universal" density functional). In this case, it is 
accepted to use the term "universal" in the sense of the independence of the explicit form of 
the external field potential [1,2], although this necessitates to consider the -y-representability 
of the inhomogeneous density. 

Here the operators T and U are the operators of the kinetic and interparticle interaction 
energies, respectively. Statement (3) is the basis of the density functional theory (DFT) 




(3) 
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widely used in various areas of physics and chemistry (see, e.g., [4,5]). However, the exact 
form of this universal functional is still unknown even for noninteracting electrons (U = 0). It 
is clear that there is one-to-one correspondence between the statements about the existence 
of the functionals v [n(r)] and F[n(r)], i.e., the existence of one of the these functionals 
predetermines the existence of another. Hence, if the functional F[n(r)] in (3) does not 
exist, the functional v [n(r)] also does not exist [3]. 

In this connection, let us consider one non-relativistic electron of mass m in the static 
external field v(r). Then the electron steady state characterized by a certain set of quantum 
numbers a, including the spin quantum number a, is completely defined by wave function 
$ a (r) which satisfies the Schrodinger equation 



where e a is the electron energy in a corresponding state. Since the electron energy is inde- 
pendent of spin, each value of e a is doubly degenerate in the spin quantum number, as well 
as other values of physical quantities, including the inhomogeneous density n a (r). It is ac- 
cepted to solve equation (4) for eigenvalues when the boundary condition 3> Q (|r| — > oo) = 
(the so-called condition at infinity [6]) is satisfied. Taking into account the possibility of 
considering the system in a finite volume V, the boundary condition for Eq. (4) in the most 
general form is written as 



where S is the surface bounding the volume V. We further take into account that the wave 
function $ Q (r) can be considered as a real function [6]. Then 



V r n a (r) = 2$ a (r)V r $ a (r), . A r n a (r) = 2$ a (r)A r $ a (r) + 2(V r $ a (r))(V r $ a (r)). (7) 



It immediately follows from (4)- (7) that the inhomogeneous density n a (r) satisfies the equa- 
tion for eigenvalues 




(4) 



$«(r -> S) = 0, 



(5) 




(6) 
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h 2 h 2 

- —A r n a (r) + rT (V r n Q (r))(V r n a (r)) + u(r)n a (r) = e a n a (r) (8) 

4m 8ron a (r) 

with the boundary conditions 

n a (r^S) = 0, V r n a (r)| r ^ s = 0. (9) 

Now let us integrate Eq. (8) over the volume occupied by the system, taking into account 
the normalization condition immediately following from (6), 



/ 



n a (r)dV = 1 (10) 
Then, from (8), we find the density functional for the energy e a 

e a [n a (r),v(r)} = F«[n Q (r)] + J v(r)n a (r)dV, (11) 



F«[n Q (r)] = F (1) [n a (r)] + F$[n a (v)}, (12) 



^kW] - / ^ F&Mr)] = £ / (V^(r))(V r n a(r) ) ^ 

Here F^[n Q (r)] is the universal density functional _F[n(r)] (3) for one particle (index (1)). 
The functional F^[n a (r)] is written in (12) in the form of two terms by two reasons. First, 
taking into account the Gauss formula and the second boundary condition in (9), the func- 
tional F ^[n Q (r)] vanishes, 

*2 p +2 p 

F (1) K*(r)] = -—J A r n a (r)dV = --|V r n a (r)« = 0. (14) 

Second, the functional explicit form F^[ra a (r)] (13) formally is exactly identical to the so- 
called Weizsaker correction to the Thomas-Fermi kinetic energy functional [7] (see [4,5] for 
more details). In contrast to the Weizsacker correction, the expression F^[n a (r)] (13) in 
this problem is an exact expression for the "universal" functional of the inhomogeneous 
density for one particle. Thus, in the case at hand, the "universal" density functional 
F^[n a (r)] exists, is found exactly, and is written as 

FW[n a {T)] = F$[n a (r)] (15) 
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By direct calculation (see, e.g., [5]), it is easy to verify that 

-s^r = -s^w A '"» (r ) + 8^iw (v '"« (r ) )(v '' l °( r ^ < 16) 

Thus, Eq. (8) for energy eigenvalues of one particle in the external field v(r) with bound- 
ary conditions (9) is a consequence of the variation equation for the energy e^[n^\r),v(r)] 
as the inhomogeneous density functional n^(r) of one particle in the specified external field 
v(r), 

5e (1) [n (1) (r)] = (17) 

Indeed, using normalization condition (10) and the Legendre transform, from (17), we find 

<5e (1, [» (1, (r)l 

fa»)(^r = const - (18) 

To determine the constant in Eq. (18), we take into account that, according to (11)-(16), 

5F«[n«(r)] . . m 

Thus, the variation equation (17) is equivalent to (8) and is completely identical to the 
corresponding equation of the particle system as the wave function functional in quantum 
mechanics (see, e.g., [6]). Hence, in the case of one particle under consideration, there exists 
the density functional v [n(r)] for the external field potential, determined within a constant 
factor, 

„[„«(,)] + const = = ^^Kn^r) - _JL_ F( v r „<-»(r))(V r n<')( r)5 20) 

Thus, as noted above, the existence of the functional F^[n^(r)] (15) predetermines the 
existence of the functional w[n(r)] (20), and vice versa. 

Taking into account that the Hohenberg-Kohn lemma proof, as well as statements (1) 
and (3), are by no means independent of a particular value of the number of particles 
iV in the system under study, the assumption on the existence of the "universal" density 
functional makes it possible to extend the results obtained to the case of an arbitrary number 
of noninteracting electrons. Then the explicit form of the functional for any number of 
particles is obtained using the following replacements: 

n (1) (r) -»n(r), e (1) ->■ E . (21) 
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In this case, the normalization condition (10) is replaced by the corresponding condition in 
(2), and its is also accepted that, according to (3), (11)-(13), that the quantity (^ol^l^o) 
for the system of several noninteracting particles is written as 

<*o|r|* > = F<°VM] =-ff A, n (r)dV + ff ggWjggWW . (22) 

Am J 8m J n(r) 

To ascertain, are relations (21) and (22) valid, we consider a system of N noninteracting 
electrons in the external field v(r). To account for the identity of electrons, it is the most 
convenient to perform such a consideration in the secondary quantization formalism (see [4- 
6] for more details). We note that such consideration is fully equivalent to the use of Slater 
determinants to describe the wave function of the system of noninteracting electrons and 
to implement the Young scheme [5,6]. Then, any state of the system of N noninteracting 
identical electrons is characterized by a set of the so-called "occupied" single-particle states 
ai,..., a n (see (4)) and, by virtue of the Pauli principle, 

ai ocj for % ^ j. (23) 

Then the energy and inhomogeneous density n^(r) in a corresponding state are given 
by [5,6] 

E {0) (ai,..., a N ) = ^e Qi , n (0) (r, a u . . . , a N ) = n Qi (r). (24) 

ai an 

In this case, e ai = e ai [n ai (r), v(r)] (see (11)-(13)), and S ai l = N. Then, it immediately 
follows from (23) and (24) that the energy of the system of iV(iV > 3) noninteracting 
identical electrons, including the ground state energy E^\ cannot be the density functional 
n (o) ( r ) j n the specified external field due to the nonlinearity of the functional F$ [n ai (r)] (13). 
It is clear that a similar statement also takes place for the universal functional (^ol^l^o)- 
This is a consequence of the fact that the external field potential v (r) at iV > 3 cannot be 
presented as the density functional n(r) (see (20)). 

In the case of two noninteracting ground-state electrons, the DFT statements remain valid 
due to the double degeneracy in the spin quantum number (n(r) = 2n^(r), see, for example, 
[8]). We also note that these results are directly associated with the Pauli principle which is 
related only to fermions. In the case of ground-state bosons which are " accumulated" at one 
lowest energy level (Bose condensation), the DFT statements remain valid at an arbitrary 
number of noninteracting bosons. 
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Thus, for more than two noninteracting fermions in the inhomogeneous ground state, the 
external field potential v(r) is not the density functional n(r), i.e., v(r) ^ v[n(r)]. 

Let us show now that the obtained results lead straightforward to validation of the anal- 
ogous conclusion for the system of interacted electrons. The Hamiltonian of the system of 
inhomogeneous electrons can be written as 



H X = T + 



J drv(r)n(r) + XU. (25) 



where n(r) is the operator of the electron density, (^o I ™( r ) I ^o) = n ( r ) an d U is the 
operator of the electron interaction energy. Then the energy E (X) of the ground state 
^o(X) satisfies to the equality (see, e.g., [9]) 

= ( * o(A) 1 d ~w 1 * o(A)) = ( * o(A) 1 u 1 * o(A)) - (26) 

From (26) straightforward follows the explicit relation for the ground state energy of the 
system of interacting electrons 

*-^"-jff <l °»- (27) 

where is the ground state energy of the inhomogeneous electron gas without interaction, 
which has been calculated above, (U)\ is the average potential energy of inhomogeneous 
electron system with the Hamiltonian (25). Take further into account that if the density 
functional i>[n(r)] exists, the average kinetic energy, as well as average potential energy, of 
the inhomogeneous electron system is the universal density functional (singly) [10]. It is 
obvious, that presence of the parameter A in the Hamiltonian (25) has no any influence 
on this statement as well as the integration in Eq. (27). Then from (27) directly follows 
that for the existence of the universal density functional F[n(r)] (see (3)) it is necessary the 
existence of the universal density functional for the inhomogeneous electron system without 
interaction. This is impossible, as it was shown above, for more then two Fermi particles. 
In a different way the impossibility of existence of the density functional of inhomogeneous 
electron system follows from the fact that the terms in the row of the perturbation theory 
for energy, which contain interaction potential, have another nature than the terms without 
interaction. Therefore, they cannot compensate the "non- universality" of the kinetic energy 
of non-interacting electrons. 

As a result, we come to the conclusion that the Hohenberg-Kohn lemma [1,2] cannot be a 
justification of the existence of the "universal" density functional as an exact statement or a 
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theorem. At the same time, in any approximations (e.g., in the limit of weak inhomogeneity 
of the external field or in the quasi-classical limit for the electron gas), the "universal" 
functional can exist. In this connection, let us pay attention to the following circumstance. 
The main approximation for density functional construction is the so-called Local Density 
Approximation (LDA) (for detail, see [4,5,10]). The basis of LDA is the dependence of the 
energy of homogeneous electron gas on average density n which equals to the ratio of the 
full number N of electrons to the volume V, n — N/V. Consideration of the homogeneous 
electron gas, it one's turn, based on use of the thermodynamic limit transition N — > oo, V — > 
oo, N/V ->n^0 (see, e.g., [9]). This means that the model of homogeneous electron system 
cannot be used as the initial approach for consideration of a finite quantity of electrons in 
an external filed (in particular, in the case of electrons in the field of one or several nuclei, 
when the conditions of the thermodynamic limit transition are not valid even for a large 
value Z» 1, where Z is the nuclear charge). 
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